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This survey gives a brief overview of theoreti- 
cally and practically relevant algorithms to compute 
geodesic paths and distances on three-dimensional 
surfaces. The survey focuses on polyhedral three- 
dimensional surfaces. 



1 Introduction 

Finding shortest paths and shortest distances be- 
tween points on a surface S in three-dimensional 
space is a well-studied problem in differential ge- 
ometry and computational geometry. The shortest 
path between two points on S is denoted a geodesic 
path on the surface and the shortest distance between 
two points on S is denoted a geodesic distance. In 
this survey, we consider the case where a discrete 
surface representation of S is given. Namely, S is 
represented as a polyhedron P in three-dimensional 
space. Since discrete surfaces cannot be differen- 
tiated, methods from differential geometry to com- 
pute geodesic paths and distances cannot be applied 
in this case. However, algorithms from differential 
geometry can be discretized and extended. Further- 
more, the discrete surface can be viewed as a graph in 
three-dimensional space. Therefore, methods from 
graph theory and computational geometry have been 
applied to find geodesic paths and distances on poly- 
hedral surfaces. 

The general problem of computing a shortest path 
between polyhedral obstacles in 3D is shown to be 
NP hard by Canny and Reif using reduction from 3- 



SAT [1 .]. Computing a geodesic path on a polyhe- 
dral surface is an easier problem and it is solvable in 
polynomial time. 

Computing geodesic paths and distances on poly- 
hedral surfaces is applied in various areas such as 
robotics, geographic information systems (GIS), cir- 
cuit design, and computer graphics. For example, 
geodesic path problems can be applied to finding the 
most efficient path a robotic arm can trace without 
hitting obstacles, analyzing water flow, studying traf- 
fic control, texture mapping and morphing, and face 
recognition. A survey related to geodesic paths in 
two- and higher-dimensional spaces can be found in 
the Handbook of Computational Geometry [31]. 

Note that the geodesic distance between any two 
points on P can be easily determined if the geodesic 
path is known by measuring the (weighted) length of 
the geodesic path. Hence, we will only consider the 
problem of computing geodesic paths on P. 

Problems on finding geodesic paths and distances 
depending on the number of source and destination 
points have been studied. The three most commonly 
studied problems are (a) finding the geodesic path 
from one source vertex s G P to one destination ver- 
tex d G P, (b) finding the geodesic paths from one 
source vertex s G P to all destination vertices in P, 
or equivalently, finding the geodesic paths from all 
source vertices in P to one destination vertex d G P, 
known as single source shortest path (SSSP) prob- 
lem, and (c) finding the geodesic paths between all 
pairs of vertices in P, known as all-pairs shortest 
path (APSP) problem. 
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The algorithms reviewed in this survey are com- 
pared by means of the following five categories: 

• Accuracy of the computed geodesic path. 

• Cost metric used to compute the geodesic path. 
The cost metric can be the Euclidean distance 
or a weight function (for example when going 
uphill is more costly than going downhill). 

• Space complexity of the algorithm. 

• Time complexity of the algorithm. 

• Applicability of the algorithm by surveying if 
the algorithm has been implemented and tested 
in practice. 

Approximation algorithms are compared accord- 
ing to their approximation ratio (or approximation 
factor) k. An algorithm that finds approximations to 
a geodesic path with approximation ratio k returns 
a path of length at most k times the exact geodesic 
path. 

To solve the problem of computing geodesic 
paths on discrete surfaces, two different general ap- 
proaches can be used. First, the polyhedral surface 
can be viewed as a graph and algorithms to com- 
pute shortest paths on graphs can be extended to 
find geodesic paths on polyhedral surfaces. Algo- 
rithms following this approach are reviewed in Sec- 
tion 2. Second, the polyhedral surface can be viewed 
as a discretized differentiable surface and algorithms 
from differential geometry can be extended to find 
geodesic paths on polyhedral surfaces. Algorithms 
following this approach are reviewed in Section 3. 
At the end of each section, open problems related to 
the section are summarized. 

2 Graph-Based Algorithms 

This section reviews algorithms to compute geodesic 
shortest paths that can be viewed as extensions of 



graph theoretic algorithms. To obtain a good under- 
standing of the reviewed algorithms, we first review 
some well-known graph theoretic algorithms. 

Dijkstra proposed an algorithm to solve the SSSP 
problem on a directed weighted graph G(V, E) with 
n vertices, m edges, and positive weights [ ]. Dijk- 
stra's algorithm proceeds by building a list of pro- 
cessed vertices for which the shortest path to the 
source point s is known. The algorithm iteratively 
decreases estimates on the shortest paths of non- 
processed vertices, which are stored in a priority 
queue. In each iteration of the algorithm, the clos- 
est unprocessed vertex from the source is extracted 
from the priority queue and processed by relaxing 
all its incident edges. The notion of relaxation under- 
lines the analogy between the length of the shortest 
path and the length of an extended tension spring. 
When the algorithm starts, the length of the short- 
est path is overestimated and can be compared to 
an extended spring. In each iteration, a shorter 
path is found, which can be compared to relaxing 
the spring. Although the original implementation 
requires 0(n 2 ) time, the running time can be de- 
creased to 0(n log n + m) by using Fibonacci heaps. 
Thorup [45] presented an 0(m) -time algorithm in 
case where each edge is assigned a positive integer 
weight. The main idea is to use a hierarchical buck- 
eting structure to avoid the bottleneck caused by sort- 
ing the vertices in increasing order from s. 

The length of a path on S depends on the em- 
ployed cost metric. Hence, the shortest or geodesic 
path on S depends on this cost metric. In Section 2.1, 
geodesic path algorithms with Euclidean cost metric 
are reviewed. Using the Euclidean cost metric im- 
plies that the Euclidean length of the path is used 
to measure the length of the path. In Section 2.2, 
geodesic path algorithms on weighted surfaces are 
reviewed. Using a weighted cost metric implies 
that different faces of S can be weighted differently. 
Clearly, any algorithm that can solve a shortest path 
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problem using a weighted cost metric can also solve 
the same problem using the Euclidean cost metric. 

2.1 Euclidean Cost Metric 

When using the Euclidean distance along a poly- 
hedron P as cost metric, shortest paths consist of 
straight line segments that cross faces of the polyhe- 
dron. An approach to compute shortest paths on P 
aims to compute a superset of all the possible edges 
of shortest paths on P and to use this information to 
compute shortest paths. Since all of the algorithms 
pursuing this strategy are mainly of theoretical in- 
terest to establish bounds on the number of possible 
edges of shortest paths on P, they are not discussed 
in detail in this survey. The algorithm's pursuing this 
approach are less efficient than the ones surveyed. 
A good overview of the algorithms finding edge se- 
quences is given by Lanthier [25, p. 30-35]. 

We first review algorithms that only operate on the 
surface of a convex polyhedron. Second, we review 
algorithms that operate on the surface of any (convex 
or non-convex) polyhedron. 

Convex Polyhedra 

This section discusses algorithms that operate on the 
surface of a convex polyhedron P with n vertices. 
Shortest paths according to the Euclidean cost metric 
are considered. 

Shark and Schorr [42] proposed an algorithm that 
computes the exact shortest path between points on 
the surface of P. The proposed algorithm is mainly 
based on three observations. First, any shortest path 
intersecting an edge e of P enters and leaves e under 
the same angle. Second, no shortest path on a convex 
polygon P can pass through a vertex p of P unless 
p is the destination or source of the shortest path. 
Third, if the sequence of edges of P intersected by 
the shortest path between s and p is known, the short- 



est path can be computed as the straight line joining 
s and p after unfolding the faces adjacent to the edge 
sequence to a plane. The authors aim to subdivide P 
with respect to a given source point s, such that the 
shortest path from s to any other point in P can be 
found efficiently. They define ridge points x of P as 
points that have the property that there exists more 
than one shortest path from s to i and prove that 
the ridge points can be represented by 0(n 2 ) straight 
line segments. The algorithm partitions the boundary 
of P into at most n connected regions called peels 
not containing any vertices or ridge points of P. The 
boundaries of peels contain only ridge points, ver- 
tices, and edges of P. The algorithm to construct 
the peels is similar to Dijkstra's graph search algo- 
rithm. The peels are then iteratively unfolded to the 
plane. The algorithm first preprocesses P by con- 
structing the peels with respect to s in 0(n 3 logn) 
time. The algorithm stores the computed peels in a 
tree called slice tree that can then be used to deter- 
mine the shortest path between an arbitrary point on 
P and s in 0(n) time. The slice tree data structure 
uses 0{n 2 ) space. 

Mount [34] improves the algorithm by Shark and 
Schorr both in terms of time and space complexity. 
The main observation by Mount is that the peels de- 
fined by Shark and Schorr can be viewed as Voronoi 
regions of a point set R containing the planar un- 
foldings of the source point s. Note that R contains 
at most n points per face of P because there are at 
most n peels intersecting a face of P. Mount ob- 
serves that the shortest path from s to any point x on 
P is at most the shortest path from x to any point 
r E R plus the distance between r and s along the 
planar unfolding of the path. This observation de- 
pends on the convexity of P and on the fact that all 
shortest paths unfold to polygonal chains consisting 
of straight line segments. Mount uses this observa- 
tion to prove that the Voronoi regions of R are iden- 
tical to the peels of P with s as source point. An 
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algorithm following the outline of Dijkstra's algo- 
rithm is used to compute the point set R and simul- 
taneously, the Voronoi regions of R. Using this ap- 
proach, P can be preprocessed with respect to s in 
0(n 2 log n) time and 0(n 2 ) space. The space re- 
quirement to store the data structure after building it 
can be reduced to 0(n log n) by storing 0(n) differ- 
ent but similar lists of size 0(n) each in an efficient 
way to avoid redundancy. Note that building the data 
structure still requires 0(n 2 ) space. The query time 
to compute a shortest path from an arbitrary point 
p € P to s is reduced to 0(k + log n), where k is the 
number of faces of P intersected by the shortest path 
by using an output sensitive point location data struc- 
ture. Mount [35] reduced the space requirement to 
build the data structure storing the Voronoi diagram 
to 0(n log n) by building a hierarchical structure on 
the intersections between edges of P and geodesic 
paths starting from s. The data structure stores for 
each edge eofPa tree whose leaves contain the in- 
tersections between e and geodesic paths crossing e 
in order. Common sub-trees of different edges are 
shared to reduce the space complexity. 

To avoid the high time complexity of finding 
geodesic paths, Hershberger and Suri [18] propose 
an algorithm that finds an approximate shortest path 
between two points on the surface of P. The algo- 
rithm takes only 0(n) time and has an approxima- 
tion factor of 2. The main idea of the algorithm is 
to extend the notion of bounding boxes to a gen- 
eral simplified representation of P and to compute 
the shortest path between two points on this sim- 
plified shape. To compute a shortest path between 
s and t, the faces containing the source and desti- 
nation points are extended into planes, and at most 
0{n) different planes are added to the two planes 
to obtain a wedge. The shortest paths between the 
two points are computed on each of these wedges in 
O(l) time as a wedge has constant description size. 
The shortest path that was found is used to approx- 



imate the shortest path on P. The algorithm can be 
extended to approximately solve the SSSP problem 
in 0(n log n) time. That is, starting from one source 
point, the algorithm computes approximations with 
approximation ratio 2 to all other points on P. 

Har-Peled et al. [ 1 7] extend the algorithm by Her- 
shberger and Suri to obtain an approximation ratio 
of (1 + e) for < e < 1. The algorithm is based on 
the approximation scheme by Dudley [14] that ap- 
proximates the minimum number of sets required to 
approximate every set as e-approximation. The al- 
gorithm by Har-Peled et al. proceeds by expanding 
P by a factor related to e and to the approximation 
obtained by Hershberger and Suri's algorithm. De- 
note the expanded polygon by P'. The shortest path 
between two vertices on P is approximated on a grid 
lattice between the boundaries of P and P'. Since 
P is convex and since the path is not in the inte- 
rior of P, the length of the path cannot be shorter 
than the true shortest path. The path obtained by this 
method can be projected to P while ensuring that the 
length of the path does not grow. The running time 
of the algorithm is 0(n min(-fs-, log n) + log ^) 
and hence depends both on n and e. As the algo- 
rithm by Hershberger and Suri, this algorithm can 
be extended to approximately solve the SSSP prob- 
lem. The running time of the extended algorithm is 
0(-o-(logra + logi)). Although the theory used 
by Har-Peled et al. is rather technical, the algorithm 
itself is simple. Agarwal et al. [2] improved the 
running time of the algorithm to approximate one 
shortest path by an approximation ratio of (1 + e) 
by Har-Peled et al. to 0(nlog^ + 4?). This im- 
proves the running time of the algorithm to approxi- 
mately solve the SSSP problem to 0( + -o- log n). 
Har-Paled [15] presents a further improvement of the 
running time of this algorithm. After preprocessing 
the convex polytope in 0(n) time, an (1 + e) approx- 
imation of the shortest path between two vertices is 
reported in 0(-t¥ + 40 time. This improves the 
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running time of the algorithm to approximately solve 
the SSSP problem to 0(n(l + *2fE£ + ^)). 

Recently, Schreiber and Shark [39] proposed an 
exact solution to the SSSP problem on convex poly- 
hedra in 3-dimensional space. The algorithm ex- 
tends Dijkstra's algorithm to allow continuous up- 
dates. That is, a wavefront is propagated from the 
source s along the boundary of P and the wavefront 
is updated at events that change the topology of the 
wavefront. Note that a similar technique of continu- 
ous Dijkstra updates was used in [? ]. The general 
idea of the continuous Dijkstra technique was for- 
mally described by Mitchell et al. [32] and is re- 
viewed later in this survey. An implicit represen- 
tation of the solution is computed in optimal time 
O(nlogn). The implicit representation is stored us- 
ing 0(n log n) space. Afterwards, the shortest path 
from the source to any point on P can be reported in 
0(log n + k) time, where k is the number of faces of 
P crossed by the path. 

Schreiber [ ] extends the previous approach by 
Schreiber and Sharir [39] to so-called realistic poly- 
hedra. Realistic polyhedra are defined as three 
classes of non-convex polyhedra. The first class of 
realistic polyhedra have a boundary that forms a ter- 
rain whose maximal facet slope is bounded by a con- 
stant. The second class of realistic polyhedra has 
the property that each axis parallel square with edge 
length I that has distance at least / from P is inter- 
sected by at most a constant number of faces of P. 
The third class of realistic polyhedra has the prop- 
erty that for each edge e of P of length |e|, there are 
at most a constant number of faces within shortest 
path distance 0(|e|). 

Agarwal et al. [1] propose an algorithm to com- 
pute a (1 + e) approximation of the shortest path be- 
tween two vertices that uses 0(-^) time and O(^) 
space. The approach proceeds by constructing a 
graph, computing the shortest path on this graph, and 



projecting the computed graph onto the surface of P. 
Agarwal et al. implemented and tested this algorithm 
and the algorithm by Hershberger and Suri [18] for 
artificial data sets with up to almost 100000 faces. 

The following problems related to computing Eu- 
clidean shortest paths on the surface of a convex 
polyhedron remain unsolved: 

• Can SSSP problem on convex polyhedra 
be solved in O(nlogn) time using 0(n) 
space [39]? 

• Can an efficient trade off between the query 
time and the space complexity be estab- 
lished [39]? 

General Polyhedra 

This section discusses algorithms that operate on the 
surface of a polyhedron P with combinatorial com- 
plexity n. Note that P need not be convex. Short- 
est paths according to the Euclidean cost metric are 
considered. The main problem that occurs when al- 
lowing non-convex polyhedra is that geodesic paths 
from s to t on P may pass through a vertex p of P. 

O'Rourke et al. [37] extend the algorithms by 
Sharir and Schorr [42] and Mount [ ] to obtain 
the first algorithm that finds the exact geodesic path 
between two vertices of an arbitrary polyhedron in 
polynomial time. Both the source and the destination 
point of P are considered to be vertices of P. The al- 
gorithm considers the problem in two steps. First, the 
straight-line distances between all pairs of vertices of 
P are found. This is achieved by extending the tech- 
nique to compute peels in [42]. Second, the short- 
est distance between the source and the destination 
vertex is found on the graph induced by the vertices 
of P. The algorithm takes 0(n 5 ) time to compute 
one shortest path on P. Since the complexity of the 
running time is high, the algorithm is irrelevant for 
practical purposes and has not been implemented. 
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Mitchell et al. [32] formalize the technique called 
continuous Dijkstra previously used in [34] to find 
shortest paths from a source point s on the surface of 
a convex polyhedron. The algorithm traverses the 
graph induced by P similarly to the graph explo- 
ration of a graph G in Dijkstra's algorithm. Edges of 
P behave like nodes in G. Since the distance from 
s on P to an edge e is not unique, e is labeled by 
a function describing the distance from s to e. The 
algorithm keeps track of a subdivision of e with the 
property that for two points p and q in the same re- 
gion of e, the shortest paths from s to p and from s 
to q pass through the same sequence of vertices and 
edges of P. Mitchell et al. observe that these sub- 
divisions of e resemble the peels used in [42]. How- 
ever, special care needs to be taken when comput- 
ing this subdivision, since geodesic paths emanating 
from s can pass through a vertex p of P. In this 
case, p is treated as a pseudo-source. The pseudo- 
source p is labeled by the geodesic distance from s 
to p. For any point x of P, the geodesic distance 
is the minimum of the shortest distance from s to x 
not passing through a vertex of P and the geodesic 
distance from the nearest pseudo-source of P to x 
plus the label of the pseudo-source. This observa- 
tion allows to compute the subdivision of s and to 
store for each region of the subdivision the distance 
to the nearest pseudo-source. For a given source ver- 
tex s, the algorithm computes a subdivision of P in 
0(n 2 log n) time and 0{n 2 ) space. Once the sub- 
division has been computed, the distance from s to 
any other point on P can be computed in O(logn) 
time. Reporting the shortest path between s and any 
other point on P takes 0(k + logra) time, where k 
is the number of faces of P crossed by the shortest 
path. If the algorithm is initialized with more than 
one source point, the subdivision obtained after the 
continuous Dijkstra algorithm ended represents the 
Voronoi diagram of the source points. Mitchell et 
al.'s algorithm is of theoretical interest, since the con- 



tinuous Dijkstra technique can be used with differ- 
ent update schemes to obtain new algorithms, as we 
saw for convex polyhedra [39]. Although the con- 
tribution by Mitchell et al. is rather technical, the 
algorithm is of practical interest as well. Recently, 
Surazhkhy et al. [44] implemented and tested the 
algorithm on data sets obtained using a laser-range 
scanner. Although the worst-case running time of 
the algorithm is 0(n 2 logn), Surazhsky et al. found 
the algorithm's average running time in their experi- 
ments to be much lower and suitable for objects with 
hundreds of thousands of triangles. The exact algo- 
rithm by Mitchell et al. is then modified to obtain an 
algorithm that solves the SSSP problem with approx- 
imation ratio (1 + e). Surazhsky et al. derive from 
their experiments that an average running time of 
O(ralogn) can be expected in practice for bounded 
approximations from one source point to all the other 
points of the mesh. 

Chen and Han [11] developed an algorithm to 
compute geodesic distances from a source point s 
on a non-convex polyhedron that does not use the 
continuous Dijkstra technique. The algorithm con- 
structs a tree called sequence tree that can be viewed 
as an extension of the dual graph of the tree con- 
taining ridge points used by Sharir and Schorr [42] 
to non-convex polyhedra. In the case of a convex 
polyhedron, the sequence tree T contains nodes con- 
sisting of an edge e of P, the image of s in the lo- 
cal coordinate system of the face incident to e, and 
the projection of the image onto e. Chen and Han 
prove that T has a linear number of nodes, contains 
all of the shortest paths, and can be built in 0{n 2 ) 
time. For non-convex polyhedra, T contains addi- 
tional leaves representing pseudo-sources of P (as 
defined by Mitchell et al. [32]) and the distances of 
pseudo-sources from s. This augmentation of T adds 
at most 0(n) nodes. Hence, the algorithm builds a 
sequence tree in 0(n 2 ) time and 0(n) space. After 
T was computed, the geodesic distance between s 
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and any point in P can be reported in 0(log n) time. 
The geodesic path can be reported in 0(logn + k) 
time, where k is the number of faces of P crossed by 
the path. Kavena and O'Rourke [20] implemented 
and tested the algorithm on synthetic data. The im- 
plementation confirms the quadratic time complexity 
and the linear space complexity in practice. The ex- 
periments show that roundoff errors are not a serious 
problem for this algorithm. Kavena and O'Rourke 
found the space complexity to be the bottleneck of 
the algorithm. Data sets with tens of thousands of 
points were used to test the algorithm. 

Har-Peled [16] extended his previous ap- 
proach to compute (1 + e) approximations of 
geodesic paths on realistic polyhedra [15] to 
work in case of general polyhedra. Given a 
source point s, the algorithm computes a sub- 
division of P of size 0(-log(-)) in time 
O (n 2 logn + j log (-) log (f))- In the special 
case of convex polyhedra, the preprocessing time 

becomes O ((f) 3 log (±) + ^ log (±) logn). 

After this preprocessing step, a (1 + e) approx- 
imation of the shortest path between s and any 
point p on P can be reported in O (log (f )) 
time. This implies that a (1 + e) approxima- 
tion to the SSSP problem can be obtained in 
<3(n 2 logn + f log (i) log (f) +nlog(f)) time. 

In 1997, Agarwal and Varadarajan [3] proposed 
an algorithm that answers the question of whether it 
is possible to compute an approximate shortest path 
between two points on polyhedra in sub-quadratic 
time. The proposed algorithm only works for poly- 
hedra of genus zero. Two algorithms using the 
same general technique were proposed. The first 
algorithm computes an approximation to the short- 
est path with approximation ratio 7(1 + e),e > 
in 0(n3 logs n) time. The second algorithm takes 
only 0(ns logs n) time, but the approximation ratio 
increases to 15(1 + e), e > 0. Note that the run- 



ning times of both algorithms is independent of the 
choice of e. The main idea of the algorithm is to par- 
tition the boundary of the simple polyhedron P into 
patches of faces of P. A graph Gi is constructed 
on the boundary of each patch Pi. The graphs Gi 
are merged into one graph G and the geodesic paths 
on P are approximated by the solution of Dijkstra's 
algorithm on G. Although this is the first paper to 
break the quadratic time complexity, the contribution 
is mainly of theoretic interest because the algorithm 
is involved. Hence, the algorithm has not been im- 
plemented. 

Kapoor [21] presents an algorithm that solves the 
problem of computing the exact shortest path be- 
tween a pair of points on P in sub-quadratic time. 
The algorithm follows the continuous Dijkstra tech- 
nique by Mitchell et al. [32] and propagates a wave- 
front over the surface of P stalling from a source 
point s. The algorithm maintains the wavefront 
as a collection of circular arcs with centers at s 
and pseudo-sources of P. Furthermore, the algo- 
rithm maintains all of the edges of P that have not 
yet been reached by the wavefront. The algorithm 
takes 0(n log 2 n) time and 0(n) space. According 
to O'Rourke [36], the details of the algorithm are 
"formidable". It is therefore not surprising that the 
algorithm contains some flaws [39]. 

Kanai and Suzuki [19] propose an iterative ap- 
proximation algorithm to compute the shortest path 
between pairs of points on P. The algorithm is 
based on Dijkstra's algorithm and iteratively refines 
the mesh in regions where the path can pass. The re- 
finement proceeds by placing Steiner points on edges 
of P and to repeat Dijkstra's algorithm on the aug- 
mented graph. The user gives two thresholds re- 
lated to the accuracy of the approximation. The first 
threshold defines the number of times the algorithm 
iterates. The second threshold is related to the num- 
ber of Steiner points placed on an edge of P. The au- 
thors implement the algorithm and compare it to an 
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implementation of Chen and Han's algorithm. They 
find their algorithm to outperform Chen and Han's 
algorithm both in terms of time and space complex- 
ity. 

The following problems related to computing Eu- 
clidean shortest paths on the surface of a possibly 
non-convex polyhedron remain unsolved: 

• Can the exact shortest path between a pair of 
vertices on P be computed in 0(n log n) time 
using 0(n) space? 

• Can the SSSP problem be solved in 0(n log n) 
time and 0(n) space? 

2.2 Weighted Cost Metric 

This section discusses algorithms that operate on the 
surface of a possibly non-convex polyhedron P with 
combinatorial complexity n in 3-dimensional space. 
Unlike in Section 2.1, the length of the shortest path 
is not simply measured by its Euclidean length. In- 
stead, a weight w-i is associated with each face fa 
of P. The length of a path crossing /j is its Eu- 
clidean length multiplied by The weights can 
be used to model the difficulty of the path. For ex- 
ample, it is harder to walk on an uneven terrain than 
on an asphalt road. A good overview of algorithms 
related to weighted shortest paths can be found in 
Lanthier [25]. 

Mitchell and Papadimitriou [33] present an algo- 
rithm to compute the shortest path distance between 
two arbitrary points in a planar subdivision with n 
edges. They note that shortest paths obey Snell's 
Law of refraction at edges of the subdivision. The al- 
gorithm is based on this observation and the continu- 
ous Dijkstra technique formalized in [3 ]. Therefore, 
the authors note that the algorithm can be extended 
to compute weighted shortest paths on the surface of 
possibly non-convex polyhedra. The algorithm finds 



an approximation of the shortest path with approx- 
imation ratio (1 + e) using 0(n 8 log(nJV^)) time 
and 0(n 4 ) space, where N is the largest integer co- 
ordinate of any vertex in the subdivision and ^ is 
the ratio between the maximum and the minimum 
weight. To our knowledge, this algorithm has not 
been implemented. This is not surprising, since the 
high time complexity makes the algorithm unsuitable 
for practical purposes. 

Lanthier et al. [26, 27] developed an approach to 
construct a graph that can be searched to obtain an 
approximate shortest path on P by adding Steiner 
points on each edge of P. Without loss of generality, 
the authors assume P to be triangulated. A total of 
0(n 2 ) Steiner points are added, yielding space com- 
plexity 0{n 2 ) for all of the algorithms. Four algo- 
rithms are presented. The first algorithm computes 
the shortest path between two arbitrary points on P 
by finding shortest paths in the graph containing ver- 
tices of P and the added Steiner points. The com- 
puted shortest path is at most WL longer than the 
true weighted shortest path on P, where W is the 
maximum weight and where L is the longest edge of 
P. The running time of this algorithm is 0(n 5 ). Sec- 
ond, a faster and less accurate algorithm is presented 
to compute the shortest path between two arbitrary 
points on P by computing a spanner on the graph 
containing the Steiner points and by finding a short- 
est path on the spanner. The computed shortest path 
has length at most /3(tv + WL), where 7r is the true 
weighted shortest path on P and > 1 is a con- 
stant. Third, algorithms were presented to process P 
for queries asking for shortest paths between a fixed 
source point s in P and an arbitrary point q in P. 
The query time is proportional to log n and the ac- 
curacy of the shortest path. Fourth, algorithms were 
presented to process P for queries asking for shortest 
paths between two arbitrary points in P. The query 
time is proportional to log n and the accuracy of the 
shortest path. For the query algorithms, time-space 



8 



trade off schemes are presented. The authors imple- 
mented and tested all of the algorithms on both real- 
life and synthetic data sets. Experiments showed that 
in practice, much less than n 2 Steiner points suffice 
to yield acceptable results. 

Lanthier et al. furthermore present an algorithm 
that runs in 0(n log n) time that computes a shortest 

path withing a factor of ( 1 + -=— J — ) , where Q m in 

\ Sill \Jmin J 

is the minimum interior angle of any face of P [27, 
Theorem 3.1]. 

Aleksandrov et al. [5] presented an algorithm to 
compute an approximation of a weighted geodesic 
path on arbitrary polyhedra with approximation ra- 
tio (1 + e). The algorithm is similar to [26] in that 
Steiner points are added along each edge of P. On 
each edge, m = 0(log — ) Steiner points are placed, 
where L is the length of the longest edge of P and 
r is m i n ( 2+3w/w > g) times the minimum distance 
of a vertex of P to the boundary of the union of 
its incident faces. As before, — is the ratio be- 
tween the maximum and the minimum weight. A 
graph G is computed on the Steiner points and G 
is partitioned into — sub-regions. In each sub- 
region, all shortest paths between pairs of vertices 
are computed. Furthermore, all shortest paths be- 
tween pairs of vertices on the boundaries of the 
sub-regions are computed. This computation takes 
0(nmr log r + log ^ + time. The 

graph G has complexity 0(nm 2 ), which dominates 
the space requirement of the algorithm. Note that 
since the graph is subdivided into small sub-graphs, 
the preprocessing is suitable for parallelizing the al- 
gorithm. After preprocessing, a (1 + e) approxima- 
tion of the shortest path between two arbitrary query 
points on P can be reported. To our knowledge, this 
algorithm has not been implemented. 

Aleksandrov et al. [6] extend this algorithm and 
place Steiner points on edges of P and in the inte- 
rior of faces of P. The approximation ratio remains 



(1 + e). That way, a graph G is constructed. An ex- 
tension of Dijkstra's algorithm can be run on G to 
obtain a (1 + e) approximation of shortest paths on 
P. The algorithm takes 0(f log ±(4j + log n)), for 
< e < 1 and 0(n log n) for e > 1 time to compute 
the shortest path between two arbitrary vertices on 
P. 

Sun and Reif [43] improve the algo- 
rithm by Aleksandrov et al. [6] to run in 
0{j (log - + log n) log -) time. This improvement 
is achieved by solving the SSSP problem on the 
graph enhanced by Steiner points using a new 
algorithm called Bushwhack algorithm. The Bush- 
whack algorithm is similar to Dijkstra's algorithm. 
However, the Bushwhack algorithm maintains for 
each Steiner point a small set of incident edges that 
are likely to be used in order to improve the current 
shortest path. This list of edges results in an algo- 
rithm that is faster than Dijkstra's algorithm. Sun 
and Reif implemented and tested their algorithm. 
They found that when O(Mog^) Steiner points 
are inserted per edge, the hidden constant in the 
O-notation is large. 

Aleksandrov et al. [7, 8] improve the running time 
of the algorithm to 0(A^ log j log |) by discretiz- 
ing P differently. In this algorithm, Steiner points 
are placed along the three bisectors of triangles of 
P. The practical use of the algorithms is limited, 
since a large number of Steiner points is inserted. 
Due to memory restrictions on current computers, 
this yields problems for real-life data sets. 

If the weights used for the weighted distances are 
restricted to be in the range [1, p] U {oo} , p > 1, 
Cheng et al. [12] present an algorithm to compute 
an approximation of ratio (1 + e) of the shortest path 
from s to t that runs in 0( p2 '° gp n 3 log ^) time. The 
advantage of this algorithm is that the running time 
does not depend on the geometry of P. 

Lanthier et al. [28] implemented the first paral- 
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lei algorithm to compute approximations of ratio 
(1 + e) for weighted shortest paths. As in previous 
approaches, the approach proceeds by constructing 
a graph and by computing the shortest path between 
vertices of a graph. The computation of the short- 
est paths is based on Dijkstra's algorithm and can be 
broken down into three components: preprocessing 
to find a graph G, executing Dijkstra's algorithm on 
G, and backtracking the path. The algorithm uses 
a spatial indexing structure called multidimensional 
fixed partition that achieves load balancing and re- 
duces the idle time of processors. The algorithm can 
solve the SSSP and the APSP problems. The algo- 
rithm was tested on a network of workstations, on a 
beowulf cluster, and on a symmetric multiprocessing 
architecture. The tests were performed for six geo- 
graphic data sets with up to one million triangles and 
achieved acceptable running times. 

Aleksandrov et al. [ ] preprocess P in 
O log \ log f ) time and O (j± log jj space, 
such that an approximation of ratio (1 + e) between 
any point q on P and a given source point s on P 
can be computed in O (-) time. Alternatively, P 

can be preprocessed in O ^ ~S75~ l°g 7 l°g 4 7) 

time and O ^ ^t/j" log 4 7) space, such that an 
approximation of ratio (1 + e) between any pair 
of points on P can be computed in 0{q) time, 
where g is the genus of P and where q is an input 
parameter. The algorithm is complex and has not 
been implemented to our knowledge. 

The following problems related to computing 
weighted shortest paths on the surface of a possibly 
non-convex polyhedron remain unsolved: 

• No exact algorithm for computing weighted 
shortest paths exists to our knowledge. 

• How can m Steiner points be placed on each 
face such that the best approximation accuracy 



is obtained [25]? Is this problem NP-hard? 

• What is the minimum number of Steiner points 
needed on each face to obtain a (1 + e)- 
approximation scheme? Is this problem NP- 
hard? 

3 Sample-Based Algorithms 

This section reviews algorithms for computing short- 
est paths on discretized smooth surfaces. We focus 
on the case where the discretization at hand is given 
as polyhedron. Unlike the above-mentioned algo- 
rithms, the algorithms reviewed in this section gener- 
alize algorithms from differential geometry to com- 
pute geodesic paths on smooth surfaces to operate 
on discretized surfaces. The research area concerned 
with these problems is discrete differential geometry. 
For a more extensive survey, refer to Kirsanov [24]. 

Kimmel and Kiryati [22] assume that a discretized 
surface is given in a voxel representation. That is, 
space is divided into a cubical grid and each grid 
point is labeled as located inside the surface, on the 
surface, or outside of the surface. The approach pro- 
posed by Kimmel and Kiryati has two stages. In the 
first stage, a 3D length estimator is used in a graph 
search on the graph defined by the surface voxels to 
find a global approximation of the shortest path. This 
approximation is then refined using local informa- 
tion. The refinement is done using a discrete version 
of geodesic curvature shortening flow. This way, an 
approximation of a shortest path between two grid 
points can be found. The approximation ratio is not 
shown to be bounded. However, since the underly- 
ing surface is assumed to be smooth, the approxima- 
tion is the best that can be obtained with the available 
voxel grid size. The algorithm has been implemented 
and tested thoroughly and appears to perform well in 
practice. 
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Kimmel and Sethian [23] present an approach 
called fast marching method on triangular domains 
(FMM) that solves the SSSP problem by solving the 
Eikonal equation on a triangular grid with n vertices. 
The result is based on Sethian's method to solve the 
Eikonal equation on a quadrilateral grid [40, 41]. 
The algorithm's running time is 0(n log n). The al- 
gorithm proceeds by iteratively unfolding all of the 
triangles of the triangular mesh. When unfolding tri- 
angles, Steiner points are placed along the edges of 
the triangles. This results in shortest paths that cut 
through faces of the triangulation and yields consis- 
tent results. However, the shortest paths found using 
the FMM method only approximate the true geodesic 
distances on the triangular mesh. The reason is that 
the true geodesic distance may require Steiner points 
in the interior of a triangular face. The accuracy of 
the approach depends on the quality of the underly- 
ing triangulation; namely on the longest edge and the 
widest angle in the triangular mesh. The algorithm 
requires 0(n) space. Since the algorithm is easy 
to implement and performs well in practice, several 
implementations of the algorithm exist. Yatziv et 
al. [47] improve the running time of FMM by us- 
ing an untidy priority queue. Their experimental re- 
sults show that the accuracy of the computed short- 
est paths only suffers slightly from this newly intro- 
duced inaccuracy. Kirsanov [24] introduces a novel 
update rule for FMM during the march. This update 
rule yields a higher accuracy of the resulting shortest 
paths. Bertelli et al. [9] consider solving the APSP 
problem using FMM. Their goal is to take advan- 
tage of the redundant computation in different passes 
of the SSSP algorithm to obtain a more efficient ap- 
proach than simply running FMM n times with each 
vertex as source point. Although the algorithm is 
shown to achieve higher efficiency in experiments, 
the worst case running time of the algorithm remains 
0{n 2 log n). 

Martinez et al. [29] present a way to iteratively 



improve an existing estimate of a geodesic path be- 
tween two vertices of a triangulated surface. Start- 
ing from a path computed via FMM, the path can 
be refined to yield a better approximation. Simi- 
lar to Kimmel and Kiryati [22], a discrete geodesic 
curvature flow is used to iteratively improve the ap- 
proximation. Martinez et al. show that the iterative 
scheme converges to the true geodesic path. 

Xin and Wang [46] present another iterative 
method to improve the path found by the fast march- 
ing method. The algorithm first improves the ini- 
tial fast marching method by classifying the edges of 
P into different types and by treating different edge 
types differently during the wave front propagation. 
Second, the algorithm iteratively improves the result- 
ing shortest path until the exact locally shortest path 
is found. 

Memoli and Sapiro [30] approximate the geodesic 
distances of an underlying smooth manifold using 
a cloud of sample points without aiming to recon- 
struct a polyhedron representing the surface. The al- 
gorithm is based on a previous algorithm that oper- 
ates on implicit surfaces. The algorithm proceeds by 
placing a ball around each sample point and by com- 
puting the union U of those balls. The Euclidean 
distance in U is used to approximate the geodesic 
distance on the underlying smooth manifold. The 
approximation is proven to be bounded by a constant 
if the sampling rate is sufficiently small. The sam- 
pling rate needed by the algorithm depends on the 
highest principal curvature of the underlying smooth 
surface. If the sampling is subject to noise, the bound 
on the approximation error gets worse. However, the 
decline in accuracy can be bounded by a function de- 
pending on the sampling noise. 
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Table 1: Results on Shortest Paths on a Polyhedral Surface P with n vertices. The constant e > is the 
desired accuracy of the shortest path. In the weighted case, N is the largest integer coordinate of any vertex 
in the subdivision and ^ is the ratio between the maximum and the minimum weight. The symbol m 
denotes the number of Steiner points placed along one edge. The symbol r denotes mm( 2+3 ^ w , |) times 
the minimum distance of a vertex of P to the boundary of the union of its incident faces. The constant p > 1 
is the largest weight assigned to a face of P. 
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The following problems related to sample-based 
geodesic computations remain unsolved: 

• Graph-based algorithms find globally optimal 
paths that may not be locally optimal if the 
graph is based on samples obtained from a 
smooth surface. Algorithms from differen- 
tial geometry can be discretized to find locally 
shortest paths. However, these algorithms can 
often get trapped in local insignificant minima. 
Can graph-based algorithms be combined with 
algorithms from discrete differential geometry 
to yield efficient globally convergent algorithms 
to compute a bounded approximation of the 
geodesic distance on a sample set obtained from 
a smooth surface [24]? 

• Can FMM be generalized to solve the APSP 
problem in o(n 2 log n) time (recall that a solu- 
tion in 0(n 2 logn) was suggested [9])? 

4 Summary 

To summarize this survey, Table 1 gives the reviewed 
results on shortest path problems on polyhedral sur- 
faces. 

References 

[1] Pankaj Agarwal, Sariel Har-Peled, and 
Meetesh Karia. Computing approximate short- 
est paths on convex polytopes. Algorithmica, 
33(2):227-242, 2002. 

[2] Pankaj Agarwal, Sariel Har-Peled, Micha 
Shark, and Kasturi Varadarajan. Approximat- 
ing shortest paths on a convex polytope in three 
dimensions. Journal of the ACM, 44(4): 567- 
584, 1997. 



[3] Pankaj Agarwal and Kasturi Varadarajan. Ap- 
proximating shortest paths on an nonconvex 
polyhedron. In IEEE Symposium on Foun- 
dations of Computer Science, pages 182-191, 
1997. 

[4] Lyudmil Aleksandrov, Hristo Djidjev, Hua 
Guo, Anil Maheshwari, Doron Nussbaum, and 
Jorg-Rudiger Sack. Approximate shortest path 
queries on weighted polyhedral surfaces. Math- 
ematical Foundations of Computer Science, 
4162:98-109, 2006. 

[5] Lyudmil Aleksandrov, Mark Lanthier, Anil 
Maheshwari, and Jorg-Rudiger Sack. An e- 
approximation for weighted shortest paths on 
polyhedral surfaces. In Scandinavian Work- 
shop on Algorithm Theory, pages 11-22, 1998. 

[6] Lyudmil Aleksandrov, Anil Maheshwari, and 
Jorg-Rudiger Sack. Approximation algorithms 
for geometric shortest path problems. In 
STOC '00: Proceedings of the thirty-second 
annual ACM symposium on Theory of comput- 
ing, pages 286-295, 2000. 

[7] Lyudmil Aleksandrow, Anil Maheshwari, and 
Jorg-Rudiger Sack. An improved approxima- 
tion algorithm for computing geometric short- 
est paths. Lecture Notes in Computer Science, 
2751(l):246-257, 2003. 

[8] Lyudmil Aleksandrow, Anil Maheshwari, and 
Jorg-Rudiger Sack. Determining approximate 
shortest paths on weighted polyhedral surfaces. 
Journal of the ACM, 52(l):25-53, 2005. 

[9] Luca Bertelli, Baris Sumengen, and B. S. Man- 
junath. Redundancy in all pairs fast marching 
method. In IEEE International Conference on 
Image Processing 2006 (ICIP), 2006. 



13 



[10] John Canny and John Reif. New lower bound 
techniques for robot motion planning prob- 
lems. In IEEE Conference on Foundations of 
Computer Science, pages 39-48, 1987. 

[1 1] Jindong Chen and Yijie Han. Shortest paths on 
a polyhedron; Part I: computing shortest paths. 
International Journal of Computational Geom- 
etry & Applications, 6(2): 127-144, 1996. 

[12] Siu-Wing Cheng, Hyeon-Suk Na, Antoine Vi- 
gneron, and Yajun Wang. Approximate shortest 
paths in anisotropic regions. SIAM Journal on 
Computing, 38(3):802-824, 2008. 

[13] Edsger. W. Dijkstra. A note on two problems 
in connexion with graphs. Numerische Mathe- 
matik, 1:269-271, 1959. 

[14] R. M. Dudley. Metric entropy of some classes 
of sets with differentiable boundaries. Journal 
of Approximation Theory, 10:227-236, 1974. 
Erratum in Journal of Approximation Theory 
26:192-193, 1979. 

[15] Sariel Har-Peled. Approximate shortest paths 
and geodesic diameters on convex polytopes in 
three dimensions. Discrete and Computational 
Geometry, 21:216-231, 1999. 

[16] Sariel Har-Peled. Constructing approxi- 
mate shortest path maps in three dimensions. 
SIAM Journal on Computing, 28(4): 1 182— 
1197, 1999. 

[17] Sariel Har-Peled, Micha Sharir, and Kasturi R. 
Varadarajan. Approximating shortest paths on 
a convex poly tope in three dimensions. In SCG 
'96: Proceedings of the twelfth annual sympo- 
sium on Computational geometry, pages 329- 
338, 1996. 



[18] John Hershberger and Subhash Suri. Practical 
methods for approximating shortest paths on a 
convex poly tope in R 3 . In SODA '95: Proceed- 
ings of the sixth annual ACM-SIAM symposium 
on Discrete algorithms, pages 447-456, 1995. 

[19] Takashi Kanai and Hiromasa Suzuki. Approx- 
imate shortest path on a polyhedral surface 
and its applications. Computer Aided Design, 
33(11):801-811, 2001. 

[20] Biliana Kaneva and Joseph O'Rourke. An im- 
plementation of Chen & Han's shortest paths 
algorithm. In Proceedings of the 12th Cana- 
dian Conference on Computational Geometry, 
pages 139-146, 2000. 

[21] Sanjiv Kapoor. Efficient computation of 
geodesic shortest paths. In STOC '99: Proceed- 
ings of the thirty-first annual ACM symposium 
on Theory of computing, pages 770-779, 1999. 

[22] Ron Kimmel and Nahum Kiryati. Finding 
shortest paths on surfaces by fast global ap- 
proximation and precise local refinement. In- 
ternational Journal of Pattern Recognition and 
Artificial Intelligence, 10(6):643-656, 1996. 

[23] Ron Kimmel and James Sethian. Comput- 
ing geodesic paths on manifolds. National 
Academy of Sciences of the USA, 95:8431- 
8435, 1998. 

[24] Danil Kirsanov. Minimal Discrete Curves And 
Surfaces. PhD thesis, Harvard University, Ap- 
plied Mathematics, 2001. 

[25] Mark Lanthier. Shortest path problems on poly- 
hedral surfaces. PhD thesis, Carleton Univer- 
sity, School of Computer Science, 1999. 



14 



[26] Mark Lanthier, Anil Maheshwari, and Jorg- 
Riidiger Sack. Approximating weighted short- 
est paths on polyhedral surfaces. In Proceed- 
ings of the 1 3 th Annual ACM Symposium on 
Computational Geometry, 1991. 

[27] Mark Lanthier, Anil Maheshwari, and Jorg- 
Rudiger Sack. Approximating weighted short- 
est paths on polyhedral surfaces. Algorithmica, 
30(4), 2001. 

[28] Mark Lanthier, Doron Nussbaum, and Jorg- 
Rudiger Sack. Parallel implementation of geo- 
metric shortest path algorithms. Parallel Com- 
puting, Elsevier, 29, 2003. 

[29] Dimas Martinez, Luiz Velho, and Paulo Cezar 
Carvalho. Geodesic paths on triangular 
meshes. SIBGRAPI: Brazilian Symposium 
on Computer Graphics and Image Processing, 
00:210-217, 2004. 

[30] Facundo Memoli and Guillermo Sapiro. Dis- 
tance functions and geodesies on point clouds. 
Technical Report 1902, University of Min- 
nesota, Institute of Mathematics and its Appli- 
cations, 2002. 

[31] Joseph Mitchell. Geometric shortest paths 
and network optimizations. In: Jorg-Rudiger 
Sack and Jorge Urrutia (Editors). The hand- 
book of computational geometry. Chapter 15, 
pages 633-701. Elsevier Science, 2000. 

[32] Joseph Mitchell, David Mount, and Christoph 
Papadimitriou. The discrete geodesic prob- 
lem. SIAM Journal on Computing, 16:647- 
668, 1987. 

[33] Joseph Mitchell and Christos Papadimitriou. 
The weighted region problem: finding shortest 
paths through a weighted planar subdivision. 
Journal of the ACM, 38(l):18-73, 1991. 



[34] David Mount. On finding shortest paths in con- 
vex polyhedra. Technical Report 1495, Uni- 
versity of Maryland, Baltimore, Department of 
Computer Science, 1985. 

[35] David Mount. Storing the subdivision of a 
polyhedral surface. In SCG '86: Proceedings 
of the second annual symposium on Computa- 
tional geometry, pages 150-158, 1986. 

[36] Joseph O'Rourke. Computational geometry 
column 35. SIGACTNews, 30(2):31-32, 1999. 

[37] Joseph O'Rourke, Subhash Suri, and Heather 
Booth. Shortest paths on polyhedral surfaces. 
In Symposium on Theoretical Aspects in Com- 
puter Science, 1985. 

[38] Yevgeny Schreiber. Shortest paths on realis- 
tic polyhedra. In SCG '07: Proceedings of the 
twenty-third annual symposium on Computa- 
tional geometry, pages 74-83, 2007. 

[39] Yevgeny Schreiber and Micha Shark. An 
optimal-time algorithm for shortest paths on 
a convex poly tope in three dimensions. Dis- 
crete and Computational Geometry, 39:500- 
579, 2008. 

[40] James Sethian. A fast marching level set 
method for monotonically advancing fronts. 
Proceedings of the National Academy of Sci- 
ences, 93(4): 1591-1595, 1996. 

[41] James A. Sethian. Level Set Methods and 
Fast Marching Methods: Evolving Interfaces 
in Computational Geometry, Fluid Mechanics, 
Computer Vision, and Materials Science. Cam- 
bridge University Press, 1999. 

[42] Micha Shark and Amir Schorr. On shortest 
paths in polyhedral spaces. SIAM Journal on 
Computing, 15:193-215, 1986. 



15 



[43] Zheng Sun and John H. Reif. On finding ap- 
proximate optimal paths in weighted regions. 
Journal of Algorithms, 58(1): 1-32, 2006. 

[44] Vitaly Surazhsky, Tatiana Surazhsky, Danil 
Kirsanov, Steven J. Gortler, and Hugues 
Hoppe. Fast exact and approximate geodesies 
on meshes. ACM Trans. Graph., 24(3):553- 
560, 2005. 

[45] Mikkel Thorup. Undirected single-source 
shortest paths with positive integer weights in 
linear time. Journal of the ACM, 46(3):362- 
394, 1999. 

[46] Shi-Qing Xin and Guo-Jin Wang. Efficiently 
determining a locally exact shortest path on 
polyhedral surfaces. Computer Aided Design, 
39(12): 1081-1090, 2007. 

[47] Liron Yatziv, Alberto Bartesaghi, and 
Guillermo Sapiro. 0{n) implementation 
of the fast marching algorithm. Journal of 
Computational Physics, 212(2):393-399, 
2006. 



16 



